The quantum optical Josephson interferometer 
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The interplay between coherent tunnel coupling and on-site interactions in dissipation-free bosonic 
systems has lead to many spectacular observations, ranging from the demonstration of number-phase 
uncertainty relation to quantum phase transitions. To explore the effect of dissipation and coherent 
drive on tunnel coupled interacting bosonic systems, we propose a device that is the quantum optical 
analog of a Josephson interferometer. It consists of two coherently driven linear optical cavities 
connected via a central cavity with a single-photon nonlinearity. The Josephson-like oscillations in 
the light emitted from the central cavity as a function of the phase difference between two pumping 
fields can be suppressed by increasing the strength of the nonlinear coupling. Remarkably, we 
find that in the limit of ultra-strong interactions in the center-cavity, the coupled system maps on 
to an effective Jaynes-Cummings system with a nonlinearity determined by the tunnel coupling 
strength. In the limit of a single nonlinear cavity coupled to two linear waveguides, the degree of 
photon antibunching from the nonlinear cavity provides an excellent measure of the transition to 
the nonlinear regime where Josephson oscillations are suppressed. 

PACS numbers: 



I. INTRODUCTION 

Cavity quantum electrodynamics (QED) experiments 
based on strong coupling of a single anharmonic emit- 
ter to a cavity-mode have lead to the observation of the 
photon blockade effect where photon-photon interactions 
alter the statistics of light emitted by the cavity. Experi- 
mental results showing nonlinearities at the single-photon 
level have recently been achieved in both atomioi^ and 
solid-state QED systems^^ with single cavities. Mo- 
tivated by the success of single-cavity QED experiments, 
center of attention has now shifted to the exploration of 
the rich physics promised by strongly-correlated quantum 
optical systems in multi-cavity and extended photonic 
media. Several works have recently considered this pos- 
sibility: a photonic version of the Bose-Hubbard model 
with an array of nonlinear cavities^^, the realization of 
the Tonks-Girardeau regime of interacting bosons in a 
non-linear optical fiber—, and a photonic analogue of the 
Kondo effect in a ID waveguide^. Most of the existing 
proposals have focused so far on strongly-correlated pho- 
tonic systems in the quasi-equilibrium regime for which 
dissipation is essentially negligible. In most realistic 
cavity-QED structures on the other hand, dissipation 
does play a substantial role and it is very difficult to 
fix the photon number. When photon losses cannot be 
neglected, the system reaches a stationary state given 
by the balance of dissipation and driving. We propose 
here an optical analogue of the superconducting Joseph- 
son interferometer, which we name the quantum optical 
Josephson interferometer, revealing new features due to 
the genuine non-equilibrium interplay of coherent tunnel- 
ing and on-site interactions. 

The two variants of the device we are presenting couple 
a central nonlinear cavity to two external driving lasers 



through either two side cavities (Fig. la,b) or two waveg- 
uides (Fig. Id). The three-cavity system can be gener- 
alized to an N-cavity system with the central nonlinear 
onei£ (Fig. lc), and in the limiting case of very large 
N this reduces to the single cavity coupled to two side- 
waveguides (Fig. Id) . In both cases, the coupling to the 
side cavities (or waveguides) is a consequence of photon 
tunneling. We assume the center-cavity can be tuned 
to have a sizable single-photon nonlinearity, e.g. due to 
some radiation-matter interaction effect, be it Jaynes- 
Cummings-type dynamics (with a single atom or quan- 
tum dot in the center cavity) 1 ^, giant Kerr nonlincar- 
ityii, or OD polariton interaction (e.g., with a quantum 
well embedded in the center-cavity)^. The model dis- 
cussed here is quite general and can be realized in a va- 
riety of quantum optical systems. In the following we 
show that light emitted from the center cavity is the re- 
sult of two competing effects, tunneling and interactions, 
leading to a crossover between the coherent and strongly 
correlated regimes. In the coherent regime, photons are 
delocalized over the three cavities and the emitted light 
is strongly dependent on the phase difference between 
the two pumping lasers (Josephson oscillations). In the 
strongly correlated regime, the inhibition of photon num- 
ber occupation beyond Fock states |0) and |1) in the cen- 
ter cavity reduces the quantum coherence between the 
two outer ones. The suppression of Josephson oscilla- 
tions in the emitted light is accompanied by a crossover 
from Poissonian to sub-Poissonian photon statistics. 

The main thrust of the present work is to investigate 
the crossover from coherent to correlated regimes of this 
intrinsically non-equilibrium system by detecting light 
emitted from the center cavity. We show that photon 
correlation measurements reveal features of an interact- 
ing few-body system that are not captured by more tra- 
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pcrimcnts, ranging from the observation of a quantum 
phase transition in Josephson junction arrays^ to the 
direct demonstration of number-phase uncertainty in su- 
perconducting islands^. While the operation of both 
electronic and photonic versions is based on the quantum 
mechanical conjugation between phase and number vari- 
ables, the quantum optical Josephson interferometer is 
an intrinsically open system. Given the enormous impact 
of the Josephson devices in electronics, from metrology 
to quantum information processing, we expect that the 
quantum optical Josephson interferometer might consti- 
tute the building block for a new class of quantum optical 
devices. 



II. THEORETICAL DISCUSSION 
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In the absence of losses the three cavity setup depicted 
in Fig. la,b is described by the Hamiltonian 
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Figure 1: The systems under consideration, (a) Schematic of 
the quantum optical Josephson interferometer with three cou- 
pled cavities. The relevant quantities of the model are also 
defined, (b) A possible photonic crystal-based implementa- 
tion (with calculated cavity mode profiles), where the middle 
cavity can contain a quantum dot or a quantum well in strong 
coupling with the high-Q photonic crystal cavity mode, (c) 
Schematic of an interferometer in the limit of a large number 
of coupled linear cavities in optical contact with the central 
nonlinear cavity, in which only the edges of the system are 
pumped, (d) A possible solid-state implementation of (c) em- 
ploying a photonic crystal circuit with two side-waveguides 
coupled to the central nonlinear cavity. 



ditional transport-type measurements. In fact, the de- 
vice we propose has close analogies with a phase-biased 
Cooper pair transistor. In the latter, the critical cur- 
rent can be electrostatically modulated by changing the 
gate potential on the central island connected to two su- 
perconducting reservoir s 16 ' 17 . The interplay between co- 
herent tunnel coupling and on-site interactions in these 
systems has been utilized in many ground-breaking ex- 



H = J2 ^kplpk 

k=l 



J(p\p2 + Pth + h.c.) + Up\p\p 2 P2 
^ (i^+h.c), (1) 



fc=1.3 



written in the rotating frame with respect to the fre- 
quencies of the two pumping lasers (h ~ 1). In 
the above equation, Ak = oju — arc the detun- 
ings of the coherent pump lasers whose amplitudes are 
■Si, 3 = | -Ei, 3 1 exp{i</>i ; 3}, where |£i,3| are assumed time- 
independent (cw pumping). We will characterize the 
three cavity system by analyzing the response of the 
center-cavity as a function of the phase difference </> = 
03 — 4>\ . The two couplings J and U quantify the hopping 
strength between neighboring cavities and the nonlinear 
photon coupling in the center-cavity, respectively. The 
operators pi,3 describe noninteracting bosonic fields in 
the external cavities, i.e. free cavity photons, while the 
elementary excitations in the center-cavity are interact- 
ing polaritons denoted by p2 . Both the tunneling and 
the coherent pumping act on purely photonic degrees of 
freedom, e.g. J is due to the overlap of the photonic part 
among nearest neighbors cavities. A rigorous derivation 
of the model in Eq. ^ inevitably depends on the spe- 
cific system under consideration, and has been provided 
before in the context of atomic* 7 - or solid-state^ cavity- 
QED implementations. Readily accessible schemes for 
achieving large photon-photon interactions are based on 
resonant coupling of a single emitter to a cavity mode 
(i.e. the Jaynes-Cummings model). As wc detail in 
the Supplementary Information section ("Experimental 
Feasibility"), the principal features we obtain using the 
model ((T|) are qualitatively identical to those predicted 
by the Jaynes-Cummings model-type single-photon non- 
linearity. 

The dynamics of the full model, Eq. ((T|), is effectively 
equivalent to that of two coupled bosonic fields: one is co- 
herently driven, while the other is nonlinear (sec Methods 
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Figure 2: Numerical solutions for the three cavity system, (a) Average population in the center cavity as a function of (f> for 
different values of U /J, in the weak pumping regime (l-EI/72 = 0.1): Josephson oscillations are barely modified as U is increased, 
(b) The same quantity as in (a), but in a stronger pumping regime (l-EI/72 = 0.7): Josephson oscillations are suppressed by 
increasing the interaction strength. The inset shows the curve calculated for U/J = 10. (c) Second-order correlation function 
as a function of J and U for 7/72 = 5 under weak pumping conditions, l-EI/72 = 0.1, showing the transition from Poissonian 
(red) to sub-Poissonian (blue) light statistics. The functional dependence Uth(J) is highlighted with a dashed white line. 



section). From now on we consider for simplicity the case 
of equal detunings and resonant pumping, = A = 0. 
Losses can be taken into account within the quantum 
Master equation in Born-Markov approximation for the 
system density matrix p, which is expressed in the usual 
Lindblad forrn 2 ^. The relevant Master equation for this 
model is 

dp 3 Ik 

— = i[p, H] + y ( 2 PkPpl - PkPkP - pp{f>k) ■ (2) 

k=l 

In most of the relevant regimes, the Master equation has 
to be solved numerically. A description of the approach 
used in this work is presented in the Methods section. 
In the remaining of the text, we will assume 71,3 = 7. 
With the specific experimental settings of Figs, lb and 
d in mind, inter-cavity tunnel coupling J ~ 1 meV— , 
cavity quality factor Q ~ 10 5 , i.e. 7 ~ 0.01 meV in the 
optical/near- infrared domain^, and nonlinear to dissi- 
pation rate ratio U/"f = 10 (see Supplementary Infor- 
mation) can be realistically achieved, which makes the 
following theoretical analysis experimentally relevant. 

It is instructive to first consider the case in which there 
is no interaction (U = 0), where an exact analytical so- 
lution for the steady state of Eq. © can be obtained. 
The case of equal amplitudes of the two driving lasers 
(Ei = £3 = E) and equal losses in the three cavities 
(7 = 72) captures all the essential details of the non- 
interacting case. In steady state the average number of 
photons in the central cavity (77,2) = (P2P2) is found to 
be 



64J 2 |£| 2 



(8J 2 +7 2 ) 



2^2 C0S J" 



(3) 



This is an analog of Josephson oscillations, imprinted in 
the light emitted from the center-cavity, due to the in- 
terference between the two coherent driving fields. Two 
features of this solution are to be noticed for a compar- 
ison with the more interesting U 7^ situation treated 



below. First, the size of the oscillations is maximized at 
J ~ 7/2, as a result of an interplay of dissipation and 
interference. Moreover, while (n2{4> — 0)) is suppressed 
and eventually goes to zero for J > 7, the oscillations 
keep a cosine-like behavior as a function of cj). 

In Fig. 2, we present our numerical results for exper- 
imentally accessible observables of the system when the 
interaction is switched on (U > 0). Rescalcd quantities 
J and E are defined for the effective two-cavity model as 
outlined in the Methods section. Do the Josephson os- 
cillations in (712(0)), as measured by detecting the light 
emitted from the center cavity (Fig. lb), remain intact? 
In Fig. 2a, we plot (712 (</>)) for various values of the inter- 
action at a pumping strength of |£'|/72 = 0.1. The size 
as well as the functional form of the oscillations barely 
change as U j J is varied across a wide range of values, un- 
der these weak pumping conditions. This picture changes 
dramatically when we pump the system stronger, shown 
in Fig. 2b for \E\/ r y = 0.7. Here, the average population 
in the center-cavity can be sizeable and nonlinear effects 
are more pronounced. In contrast to the weak pumping 
case (Fig. 2a), the size of the oscillations is suppressed to 
a great extent as U is increased. Besides the strong sup- 
pression of visibility, Fig. 2b shows a dramatic deviation 
from the cosine-like functional form, Eq. ^ , as Uj J is 
increased from zero (the behavior for U/J 1 is shown 
in the inset). 

Next we investigate how this cross-over is reflected in 

the photon statistics of light emitted from the center 

cavity. For this, in Fig. 2c we plot the zero-time de- 

(2) 

lay second-order correlation function g 2 (0) (See Meth- 
ods section) as a function of the scaled quantities C//72 

~ (2) 
and J/72- We find that g 2 (0) displays a sharp tran- 
sition from Poissonian to sub-Poissonian light statistics 
as the interaction strength U is increased. The thresh- 
old for anti-bunched (sub-Poissonian) light generation, 
Uth, is a function of J. For J/72 <C 1, the anti- 
bunching threshold is Uth(J) ~ 72, while for J/72 ^ 1, 
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Figure 3: Variation of the energy levels of the three-cavity 
Josephson interferometer as a function of U/J. We plot the 
first few total photon number manifolds in the range Ntot = 1 
to Ntot = 4, and J/72 = 0.05. The various Ntot manifolds 
are marked on the vertical axis and set off from each other 
by an arbitrary ljl for visibility. The energy levels in each 
manifold undergo an anticrossing at U/J ~ 1 and a cross- 
over takes place to an effective Jaynes-Cummings sequence as 
U/J -» 00. 



Uth(J) ~ 7 + 72- These two regimes are connected by a 
smooth crossover region. This peculiar behavior of the 
anti-bunching threshold is related to the effective dissipa- 
tion rates of the coupled system as the coupling strength 
J is varied. At small J, the coupling to the center cav- 
ity is pcrturbative; the nonlinearity (i.e. antibunching) 
therefore sets in when U is larger than the broadening 
of the bare center cavity polariton states i.e. 72. As J 
is increased, the coupling becomes non-perturbative and 
the relevant eigenstates of the coupled system are super- 
positions of center and outer cavity states; such dressed 
states have broadening contributions coming from both 
center and outer cavities and therefore the nonlinearity 
now has to be larger than the broadening of the dressed 

states in order for the system to exhibit antibunching. A 

( 2) 

remarkably simple expression for g^ (0) can be derived 
in the weak pumping limit (Sec Discussion in Supplemen- 
tary Information section) which captures all the regimes 
discussed: 



r 2 



r 2 + 4a 2 (J)/J 2 



(4) 



where T = 7 + 72 and a(J) = (4 J 2 + 7T)/(4J 2 + 772). 
These results are consistent with the expectation that 
strong photon nonlinearity can lead to photon-blockade^ 
in the center cavity giving rise to anti-bunching. While 
the relative strength of U with respect to tunnel coupling 
J seems to matter at small couplings J/72 (i.e. Uth(J) is 
a monotonic function of J), at larger J the relative effect 
of U saturates. 

An interesting feature that does not leave any footprint 
in these two observables considered is the nature of the 
system's effective non-linearity. For [/« J, the devia- 



tion of the system energy levels from a harmonic struc- 
ture is linearly proportional to U , and this determines the 

( 2) ~ 
main behaviour of g\ (0). For U ^> J however, the cen- 
ter cavity acts as a two-level system (only and 1 photon 
states available), which is coupled to the linear cavities 
with strength J. Thus, the coupled system maps onto an 
effective Jaynes-Cummings (JC) model where, surpris- 
ingly, the tunnel coupling strength plays the role that is 
commonly played by the atom-cavity dipole-coupling in 
the original JC modeUi In Fig. 3, we show how this 
comes about: two levels of each constant photon- number 
manifold (Ntot = n 2 + n>s) split off from the rest of the 
levels to form a JC-sequence as U is increased beyond J. 

In the limit of an infinite number of linear cavities 
coupled on either side to the center cavity (Fig. lc), 
we obtain a band of bosonic modes which mimic two 
external waveguides. The corresponding photon cre- 
ation/annihilation operators in Eq. |T]) can in this limit 
be replaced by their average (coherent state) values, 
P1.3 — * (^1,3) = —2iEi^/"f. The effective Hamiltonian 
of the system 



H ~ A2P2P2 



UpIpIp2P2 



EcspI + h- c - 



(5) 



then reduces to that of a single nonlinear cavity pumped 

by a coherent field with amplitude E e g = —2iJ[E\ + 

E 3 )/^. We choose parameters such that |-E e ff| = J, i-c. 

J acts as the effective pumping rate. The steady state 
(2) 

results for g^ (0) and (712) are shown in Figs. 4a and 

(2) 

4b, respectively. The g\ (0) displays a sharp transition 
from Poissonian to sub-Poissonian light statistics as the 
interaction strength U is increased. The threshold for 
anti-bunched (sub-Poissonian) light generation, Uth, is a 
function of J. For J/72 <C 1, the threshold is Uth ~ 72, 
in the opposite limit Uth ~ J- At small hoppings, the 
nonlinearity (i.e. anti-bunching) sets in when U is larger 
than the broadening of the bare center cavity polariton 
states i.e. 72. Much more interesting is the fact that at 
larger hoppings the threshold scales with J, in contrast 
to the case of three cavities (Fig. 2c) . The crossover from 
bunching to antibunching behavior reflects in a clear way 
the crossover from dclocalized to localized states. As J 
is increased, the relevant eigenstates of the coupled sys- 
tem are superpositions of center and outer cavity states. 
When J ^> U photon states are delocalized over the 
whole systems while in the opposite case the states are 
Fock states due to the onset of photon blockade. One 
may also understand the dependence of the crossover on 
J by relating it to the low- J limit of Fig. 2c: from the 
perspective of the center cavity, the driven waveguide is 
analogous to a driven cavity with a dissipation rate larger 
than all other energy scales. 

Due to the Heisenberg uncertainty relation, the 
crossover from bunching to antibunching behavior man- 
ifests itself also in the phase dependence of (712) shown 
in Fig. 4b. On increasing the interaction the visibility is 
strongly suppressed and furthermore there is a marked 
deviation from the simple cosine law found for U = 0. 
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Figure 4: Numerical solutions for the waveguide-coupled limit. We assume the coherent states |(pi,3)| = 1- (a) Second-order 
correlation function at zero-time delay for light emitted from the center cavity, as a function of U and J. A sharp crossover 
between Poissonian (red) and sub-Poissonian (blue) statistics is seen in the U — J plane, (b) Average population in the center 
cavity for J/72 = 2 as a function of </>, for different values of U /J. The oscillations approach the limiting case of Eq. (6) when 
U/J 2> 1. The inset is a zoom on the curve for U/J — 10, showing (712(0 = 0)) — » 0.5 and a strong deviation from simple 
cosine-like behaviour, (c) Scaling of '(0) and (ri2(4> = 0)) as a function of U/J for different values of J/72, showing a smooth 
crossover from delocalized to localized regime for U/J > 1. 



An analytical expression for the function in the inset can 
be found in the infinite-!/ limit, where Eq. (J5J) is replaced 
by a two-level system coupled with a driven cavity mode 



x 2cos 2 (0/2) + ( 7 2 /8J|£;|) 2 ' y ' 

which agrees with the numerical results. The broadening 
of the dip at <fi — ir is proportional to 7 2 /(4J|£'|) im- 
plying that the visibility goes to zero by increasing the 
field ((n 2 ((f> = 0)) -> 0.5). The behavior of (n 2 (0)) is 
again witness of the crossover from the delocalized to the 
correlated regimes. For large hopping the state of the sys- 
tem is approximately a coherent state. Phases are locked 
(there are strong fluctuations in the number operator) 
and the visibility is large. On the opposite case, due to 
photon blockade the state is close to a Fock state. Phase 
fluctuations in the central cavity suppress the global co- 
herence of the system and the visibility is suppressed. 

We note that although the functional forms are different, 

(2) 

both \n2\4> = 0)) and g 2 (0) display a cross-over that is 
dependent only on the dimcnsionless ratio U/J, as shown 
in Fig. 4c. 

III. SUMMARY 

In conclusion, we have studied the out-of-equilibrium 
interplay of tunneling and interactions in realistic quan- 
tum optical devices. For the three-cavity system, the 
inhibition of particle number fluctuations beyond Fock 
states |0) and |1) in the center-cavity reduces quantum 
coherence between the two outer ones. This is found to 
exhibit a threshold behavior as a function of the correla- 
tion energy U, which is clearly discernible from the satu- 
ration of Josephson oscillations at amplitude (712) = 0.5, 
and by a strong deviation from cosine-like behavior at 
strong pumping. Remarkably, we find that the effec- 
tive non-linearity of the system displays a cross-over 



into a JC-likc non-linearity as U/J is increased beyond 
1. Photon correlation measurements for this device re- 
veal a sharp threshold from Poissonian to sub-Poissonian 
statistics that is almost insensitive to the strength of the 
tunnel-coupling J. On the other hand, in the case of 
a linear array of cavities coupled to a non-linear center 
cavity, the anti-bunching threshold is found to depend 
strongly on the tunnel-coupling. This observation signi- 
fies that photon correlation measurements are very effec- 
tive in revealing the interplay of coherent tunneling and 
on-site interactions and may contain the key to inter- 
pret and probe possible phases of extended cavity-arrays 
which operate under non-equilibrium conditions. 

IV. METHODS 

The model in Eq. ([T]) is reformulated by introducing 
the canonically transformed bosonic operators s = (pi + 
P3VV2 and d = (pi —ps)/y/2, from which 

H s = A{sKs+plp 2 )+J(pU+s f P2)+Upiplp2p2+E^+E*s 1 

(7) 

where we defined = A and we discarded the dynamics 
of the field d, which is decoupled from p2- Rescaled quan- 
tities are defined as J = \/2J and E = <J2{E X + E 3 )/2. 
Thus, the dynamics of the full model ([1]) is equivalent to 
that of two coupled bosonic fields: s is coherently driven, 
while p-2 is nonlinear. Losses arc taken into account 
within the quantum master equation in Born-Markov ap- 
proximation for the system density matrix, Eq. ([2]), for 
the field operators s and p2 with dissipations 7 and 72, 
respectively. 

Analytical solution. An analytical solution to the 
steady state master equation can be found in the non- 
interacting limit, U — 0. The equation of motion for 
a generic operator expectation value, (A), is d(A)/dt = 
= i{[H s ,A]) ss + (C[A]} ss , where H s = {Jp 2 +E)s^ +h.c, 
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and the Liouvillian C[A] is formally written as 

C[A] = ^(2^As-s^A-AsH) + ^(2pIAp2-pIp2A~ApIp2 

(8) 

Solving for p\ and s, respectively, we get a system of two 
coupled equations, from which the steady state solution 
is \(p 2 ) ss \ = \E\/{ J[(l +772/(4J 2 )]}, and hence Eq. ©. 

Numerical solution. Extensive numerical simulations 
for the effective model (J7|) can be performed quite effi- 
ciently and in a reduced Hilbert space with respect to 
the full model. After explicitly defining the operators in 
matrix form on a Fock basis of bosonic number states, 
the steady state density matrix for any given set of pa- 
rameters can be obtained by finding the eigenvector cor- 
responding to the zero eigenvalue of the linear operator 
equation L\p)) = A|p)), where \p)) is the density operator 
mapped into vectorial form, and L is the linear matrix 
corresponding to the Liouvillian operator in the right- 
hand side of Eq. Q 23 ' 24 . Once \p ss )) is obtained from 
L\p)) ss = A ss |p)) ss with A ss = 0, we can recast it in matrix 
form and calculate any observable we are interested in. In 
particular, in this work we deal with (n 2 ) = Tr{p\p2p ss }, 
and the steady state zero-time delay second-order corre- 
lation function g^i? = 0) = Tr{p\p\p2P2Pss\ I ( n 2) 2 - To 
check convergence with the number of Fock states in the 
basis as a function of J, numerical results for U = are 
compared to Eq. (|3|). 



V. SUPPLEMENTARY INFORMATION 

Experimental feasibility. In the main text we have con- 
sidered a generic Kerr nonlinearity as the source of strong 
photon correlation in the center cavity. A possible way 
of experimentally implementing an effective hamiltonian 
of the type (1) or (5) in the main text is to couple 4- 
level atomic ensembles with microtoroid resonators, as 
described in detail in the literature 7 -. The latter certainly 
represents an interesting possibility for a practical real- 
ization of our proposal with state-of-the art atomic cavity 
QED. Our focus here will be on the scheme represented 
in Fig. lb and d (see text), in which a quantum dot (QD) 
is assumed to be deterministically coupled to the high-Q 
photonic crystal cavity mode in the middle^. It has been 
experimentally shown that such a system displays single- 
photon nonlinearities under coherent resonant pumping^. 
Photonic crystal (PC) circuits allow for a straightforward 
on-chip implementation with side-coupled cavities (as in 
Fig. lb) or waveguides (as in Fig. Id). If the hopping pa- 
rameter J is small compared to the laser intensities \E\\ 
and I-E3I, we can approximate the states of the external 
cavities with coherent fields of intensity 2\Ei^\/~f (with 
7 being the damping parameter of the external cavities). 
Thus, to first order in J/\E\^\, the dynamics of the cen- 
tral cavity can be effectively described by replacing in 
the Hamiltonian the operators p\ and fa with 2\Ei^\/~/, 
with reference to the derivation of the model in Eq. (5) 



in the text. With this choice the 3-cavity set-up can be 
effectively reduced to a Jaynes-Cummings model coupled 
to external driving fields and described by the simplified 
model 

H ~ A c ala 2 +A x a + &-+ig(ald-^-a + a2)+E e ffal+E* s a2 , 

(9) 

where 0,2 (aj) represents annihilation (creation) of cav- 
ity photons, while and <r + are Pauli lowering 
and rising operators related to the effective two-level 
system representing the QD exciton transition; g = 
H(Tre 2 f I em*V e f[) 1 / 2 is the exciton-photon coupling (ex- 
pressed in terms of the effective cavity mode volume and 
the QD oscillator strength)^, and A c = uj c — uj-l and 
A x = uj x — o->l are the cavity and exciton detuning from 
the pump frequency, respectively. The effective pumping 
strength is \E e g\ cx \Ei + S2I/7 and is responsible for the 
effects of phase detuning, (f> — 4>i ~ fa- 

To give some numbers, state-of-the art solid state QED 
with GaAs-based materials allows for Q ~ 10 5 — 10 6 , 
i.e. realistic 7c - 0.01 meV^i, j x ~ 1.3 X 10" 3 meV-22, 
g = 0.11 meV— . In the experiment, g is fixed as well 
as J, but the effective single-photon nonlinearity can 
be tuned by changing 8 = uj x — uj c . There are a num- 
ber of different techniques to deterministically tune the 
cavity mode frequency 4 ^ 27 ' 28 ^ 9 and/or the QD exciton 
resonance 3 ^. For the system excitation, the same laser 
source can be sent through a beam-splitter, one of the 
arms going directly into the PC circuit (e.g. through a 
tapered access waveguide) with phase fa , while the other 
being delayed and sent through a second tapered waveg- 
uide into the circuit with phase fa. Given that lasers 
with sub-MHz linewidth (i.e. much smaller than cav- 
ity and exciton dissipation rates) are currently available, 
we do not regard possible phase fluctuations in the two 
driving fields MS M. limiting issue for this scheme to be 
realized. Finally, it is key to this experiment that the 
pump laser frequencies be tuned to the lower polariton 
frequency of the JC spectrum, i.e. for each detuning 8 
we set wl ~ (uj x + u c )/2 — \/ g 2 + 8 2 /4^. 

We simulate this model by solving the corresponding 
master equation that incorporates realistic cavity and 
exciton dissipation rates. Figure [5] shows that the pre- 
dictions of the Jaynes-Cummings model in Eq. |9]) arc 
qualitatively similar to the ideal Kerr nonlinearity model 
of Eq. (5) in the main text. In particular, the de- 
pendence of light intensity on <f> shows suppression of 
Joscphson-like oscillations as the exciton frequency is 
tuned (from the blue-side) in resonance with the cavity 
mode. The increase of oscillations amplitude towards the 
bare-cavity limit appears together with a crossover from 
sub-Poissonian to Poissonian statistics as a function of 
8, consistent with the results of Fig. 4 in the main text. 
Here, tuning 8 is a way of effectively tuning the nonlin- 
earity in the central cavity, and thereby to experimen- 
tally observe the crossover from tunnel-coupled to strong 
correlated photon dynamics in a state-of-the art device. 
Notice that the linear regime is completely recovered for 
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Figure 5: Numerical simulation of the experimental feasibil- 
ity. We assume parameters of the model in Eq. Q <?/7c = 10, 
g/^fx = 100 (as realistically achievable, see text), and effec- 
tive pumping strength |-E c ff|/7c = 1 at = 0. Results are 
shown for light intensity and second-order correlation func- 
tion emitted from the middle cavity, i.e. (n c ) = (a\a2) 
and g^ 2 \0) = (a\a\a,2a2) / (n c ) 2 . (a) Josephson-like oscilla- 
tions are suppressed when the cavity-exciton detuning 5 — » 0. 
When 8 S> 7c, g, i.e. the exciton resonance is strongly blue- 
detuned from the cavity mode, the lower polariton is more 
and more cavity-like, hence the effective nonlinearity of the 
system is tuned through S. In the inset, a zoom on the curve 
for 5 = 0. (b) The crossover from correlated to delocalized 
regimes is shown for (n c ) and <? (0) at <j> = 0, respectively, 
as a function of 5. 
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Figure 6: Numerical simulation for the model of Eq. (5) in the 
main text with effective pumping strength |-E c ff |/7 = 1 at </> = 
0. (a) Josephson-like oscillations are suppressed on increasing 
the effective Kerr nonlinearity, In the inset, a zoom on 

the curve for C//72 = 6, to be compared to the previous figure, 
(b) Crossover from correlated to delocalized regimes for (712) 
and <? 2 2 ^(0) at = 0, respectively, as a function of /I//72. 



S/'jc ~ 100, i.e. S ~ 1 meV with the parameters given 
above, which is perfectly within reach of present exper- 
imental capabilities. In order to compare the behaviour 
of the different models, we notice that the strongly cor- 
related limit of Fig. (S = 0, curve in the inset) is 
qualitatively similar to the result obtained with a Kerr 
nonlinear model, Eq. (5) in the text, for U = (2 — V2)g, 
i.e. U/j2 — 6, as it is shown in Fig. [5^, and its inset. The 
crossover from delocalized to localized regimes as a func- 
tion of U is shown in Fig. [6}d to be compared to Fig. [SJd. 
The strict analogies between the experimental realization 
proposed and the model considered in the main text be- 
fore is maintained as long as the pumping strength of 
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|2,0> 




Figure 7: Energy level diagram and rates for the coupled 
cavity system. The various Ntot manifolds are set off from 
each other by an arbitrary uj p for visibility. The numbering 
of the states \n) are from down to up and left to right. For 
instance |1) = |1, 0). 



the Jaynes-Cummings system is low enough to discard 
higher-lying photon manifolds. 

Finally, we stress that the same model and experi- 
mental approach can be used to study the realization 
of quantum optical Josephson interferometers based on 
different technologies, such as circuit QETA Moreover, 
we point out again that the effective Kerr nonlinearity 
can be achieved, in the same system represented in Figs, 
lb and d in the text, with a QW in strong coupling to 
the cavity mode. Although there is no experimental evi- 
dence of Kerr nonlinear behaviour of 3D confined cavity 
polaritons at time of writing, it is likely that polariton 
blockade on a single quantum box can be achieved in the 
near future along the lines and numbers quoted in the 
literature^, e.g. with an alternative cavity geometry re- 
cently realized^!. Such a result would make possible a 
more direct realization of our model with a Kerr nonlin- 
earity in the solid state. 

Derivation of Eq. (4). Here we investigate the steady 
state dynamics of the Hamiltonian (7) (see Methods sec- 
tion, main text) in the weak pumping limit. We will 
present a derivation of the expressions for {ri2(4> = 0)) 

and <?2 ( T = 0) = {P2P2P2P2) / '( n 2) 2 ■ Consider the low- 
energy excitations of the Hamiltonian (7) in the basis 
\ri2,n s } where N to t = 712 + n s is the total number of 
photons. In the weak pumping limit, the drive term 
Es' + E*s causes transitions between the manifolds Ntot 
and Ntot + 1. Let us write the total time-dependent wave- 
function of the system as ^>(t) = J2 n a n\ n )- In the weak 
pumping limit, we consider the lowest three total photon 
manifolds N to t = 0, 1, 2, hence n = 0, . . . , 6. The corre- 
sponding energy level diagram and the rates are shown 
in Fig. El 

Terms neglected in 'F(i) will be of order O(^p) 3 where 
7 is the typical decay rate of the system (72, 7). It's 
important to write the equations of motion in the bare 
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basis instead of the basis of the coupled cavity states 
to get the dissipation rates correctly. The equations of 
motion arc 



It's interesting to note that in the limit J/7 ^S> 1, \a± \ ^> 
|d2 1 j showing that interference effects play an important 
role. Thus, 



= 

01 = - 

a 2 = - 



n p ~ = 16J 2 /5 2 /(772 + 4J 2 ) 



t2\2 



72 

iA + — ) 5i — i J5.2 



z'A + —J 5,2 — i Jfii — iEclq 
fi 3 = -(2zA + 2iC/ + 72)03 -z-v/2Ja 4 

04 = — ^2iA + — \ 04 — iV2J(a,3 + 05) — i-Eai 

05 = — (2iA + 7) 05 — iv2 Ja4 — ii?a2 

We have kept terms that are the same order of magnitude 
in ^ and a„ = o ra e IWL — . The steady-state solutions can 
be easily determined by additionally employing the nor- 
malization condition J2n=i \ a n\ 2 = 1 (to order 0(^-) 6 ). 
We find, for A = 0, 



AJE 



Oi = 



772 + 4 J 2 ' 



"2 



2iE l2 
772 + 4 J 2 



(10) 



(11) 



to order O(^) 4 . Solving for the two-photon manifold 
amplitudes as well, we find (4) 

2|o 3 | 2 r 2 



9?\t = 0) 



|oi| 4 T 2 + 4a 2 (J)U 2 



(12) 
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